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Abstract: Typically, the Electromechanical Impedance (EMI) technique does not use an 
analytical model for basic damage identification. However, an accurate model is necessary 
for getting more information about any damage. In this paper, an EMI model is presented 
for predicting the electromechanical impedance of a cracked beam structure quantitatively. 
A coupled system of a cracked Timoshenko beam with a pair of PZT patches bonded on 
the top and bottom surfaces has been considered, where the bonding layers are assumed as 
a Kelvin- Voigt material. The shear lag model is introduced to describe the load transfer 
between the PZT patches and the beam structure. The beam crack is simulated as a 
massless torsional spring; the dynamic equations of the coupled system are derived, which 
include the crack information and the inertial forces of both PZT patches and adhesive 
layers. According to the boundary conditions and continuity conditions, the analytical 
expression of the admittance of PZT patch is obtained. In the case study, the influences of 
crack and the inertial forces of PZT patches are analyzed. The results show that: (1) the 
inertial forces affects significantly in high frequency band; and (2) the use of appropriate 
frequency range can improve the accuracy of damage identification. 

Keywords: electromechanical impedance; structural health monitoring; PZT; Timoshenko 
beam 
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1. Introduction 

In recent years, the Electromechanical Impedance (EMI) technique has emerged as a promising 
structural health monitoring (SHM) method. It has been successfully applied to various engineering 
disciplines, including aerospace and aircraft structures [1-6] and civil structures [7-9]. In this 
technique, a piezoelectric ceramic lead-zirconate-titanate (PZT) patch is surface bonded to the 
monitored structure or embedded into a new composite construction and excited by an alternating 
voltage sweeping signal. Any physical change in the structures will result in a change of the structural 
mechanical impedance. Due to the electromechanical coupling between the PZT transducers, the EMI 
signature extracted from PZT transducer is directly related to the mechanical impedance of the 
monitored structure. Consequently, any structural physical change will induce changes in the EMI 
signature of the PZT transducer. Hence, for SHM applications, PZT EMI spectra can be compared 
with a baseline measurement during the service period of the monitored structure. Any change in the 
spectra is an indication of a change of the structural integrity, which may be caused by the presence of 
damage. 

Although the EMI-based SHM does not typically use an analytical model for basic damage 
identification, it is necessary to establish a model for more advanced features of SHM, such as damage 
evaluation and prognosis [10]. Many researchers have investigated the EMI model. As early as the 
1990s, Liang et al [11] proposed the first one-dimensional (ID) EMI-type model for a PZT-structure 
interaction system. In his model the electromechanical admittance (inverse of impedance) of PZT can 
be expressed analytically by structure impedance. Zhou et al [12] extended the ID impedance method 
to model a two-dimensional (2D) PZT element coupled to a 2D monitored structure. On the basis of 
the concept of "effective impedance", Bhalla and Soh [13,14] improved Zhou's model. Yang presented 
a generic model for predicting the electromechanical impedance of one-dimensional and 
two-dimensional PZT — structure interaction systems [15]. Based on the concept "sum impedance", 
Annamdas and Soh [16,17] proposed a three-dimensional (3D) EMI model. In general, all the models 
above ignored the adhesive layer between the PZT transducer and the monitored structure. However, 
many experimental and theoretical analysis results [18-22] have demonstrated that the mechanical 
interaction between the PZT transducer and the monitored structure occurs through interfacial shear 
stress. Hence, the adhesive layer has to be considered. Ong et al presented an EMI model which 
considered the shear lag effect of the bond layer [23]. Suresh Bhalla et al [24] incorporated the shear 
lag effect into the existing ID and 2D EMI and obtained an improved model for them. 
Yan et al [25,26] considered the adhesive layer, and presented an EMI model for Timoshenko beams 
and Mindlin-Herrmann rods. Although the factor of the adhesive layer was considered, the inertial 
terms of PZT and adhesive layer produced by the motion along with the monitored structure have not 
been taken into account in the above models. Pietrzakowski [22] noticed the inertial terms of PZT 
patches and studied the influence of bonding layer on the beam response, but his method is not suitable 
for high frequency EMI techniques due to the use of Euler beams. 

For some simple pristine structure, the analytical EMI model can be obtained. However, when 
damages are induced in structures resulting in the possibly inhomogeneity of material properties, it is 
difficult to derive an analytical formulation. In order to quantitatively identify structural damages, 
numerical methods and approximate approaches have been adopted. Naidu and Soh [27] and Tseng 
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and Wang [28] obtained the relationships between the EMI signatures and the structural changes by 
the Finite Element Method (FEM). However, for the purpose of predicting the response accurately in a 
high frequency range, a very small size element is needed. Therefore the FEM method is 
time-consuming. Yang [29] and Xu [30] applied the p-Ritz method to establish EMI models for health 
monitoring of beams, plates and cylindrical shell structures with various boundaries. Their models 
could calculate the EMI response effectively and accurately below 50 kHz, however, the accuracy of 
prediction was reduced with the increased frequency. 

In order to quantitatively identify the structural damage in a more efficient way, a more accurate 
EMI model is needed. This study proposes an electromechanical impedance model for health 
monitoring of beam structures. Different to the existing EMI model, the proposed model not only 
contains physical parameters of damage and the shear lag effect of the bonding layer, but also includes 
the transversal inertial forces of PZT patches and adhesive induced by the transversal motion of 
beam structures. 

In this paper, a pair of PZT patches bonded symmetrically onto the top and bottom surfaces of a 
rectangular beam with a crack was activated out of phase to create a pure bending excitation. Pure 
extensions in PZT patches were assumed. A shear lag model was applied to describe the behavior of 
bonding layer which is assumed to be a Kelvin- Voigt material. The PZT-adhesive layer-cracked beam 
coupled structural system was considered. The coupled system is divided into four sections due to the 
crack cross-section and the location of PZT patches. The crack is simulated by a massless torsional 
spring. Taking into account the inertial forces of PZT patches and bonding layer, the boundary 
problem is formulated by the dynamic equations of the coupled system. According to the boundary 
conditions and continuity conditions, the solution of the coupled system can be obtained, and then the 
analytical relationship between PZT admittance (inverse of impedance) and the damage parameters 
such as location and depth is derived. Finally, numerical results are presented and discussed to validate 
the proposed theoretical analysis. 

2. Formulation of the Coupled Structure and Analysis 

As shown in Figure 1, two identical PZT patches are bonded symmetrically onto the top and bottom 
surfaces of the Timoshenko beam. Two out-of-phase alternating electrical fields are applied to the PZT 
patches. 

Figure 1. Beam with a pair PZT patches bonded on its surface. 
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The beam was activated to create a pure bending vibration. A crack with a depth of h c is located at 
jci. For a better analysis, the beam is divided into four parts due to the crack cross-section and the 
location of PZT patch (Figure 1). The dynamic behavior of each part is governed by different 
equations. According to boundary conditions and continuity conditions, the dynamic equations of 
PZT-structure interaction system can be solved. The analysis method used here is similar to the 
method used by Pietrzakowski [22], but different from that work, the Timoshenko beam theory and the 
inertial of bonding layers are taken into account in this paper. 

2.1. Dynamic Equations 

The governing equations for the section with PZT patch ( x 2 < x < x 3 ) are derived considering an 
infinitesimal element shown in Figure 2. The longitudinal motion of the PZT patches is defined by 
introducing inertial forces. The governing equations can be expressed as: 

d 2 u n d 2 u n t 

O U v U X 

where 2? , p v , h p , denote Young's modulus, density and thickness of the two PZT patches, respectively; 
u p and u' v are longitudinal displacement of the upper PZT and the lower PZT, respectively; 

The shear stresses r and r' transmitted by the piezoelements are determined by the following 
stress-strain relation: 

(u —u ) 

r a V p s / 

L — 

K 

where u s and u s are the upper and lower beam surface longitudinal displacement, respectively; 
/z a is the bonding layer thickness; /z b is the beam thickness; and w is the transverse displacement of 
beam; G* is the linear function of differential operator which for the Kelvin- Voigt model of bonding 
layer material can be expressed as: 

G>G a {l + C b ^) (3) 

where G a is Kirchhoff s modulus; Cb denotes the retardation time. 

Because the two PZT patches are identical and are actived by two out-of-phase alternating electrical 
fields with same amplitude, so, the following relationship can be obtained: 

/ = r 
u = u 

P P (4) 
h b dw 

s s 2 dx 
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Figure 2. Infinitesimal beam element with piezoelements and bonding layers. 

V////MM F + dF 




Considering the dynamic coupling between the piezoelement and the beam, and including the 
mass of piezoelement and adhesive in the inertial force of intensity qt (shown in Figure 2), the 
transverse motion of Timoshenko beam is described. The governing equations of beam are 
expressed as follows: 



dM , , ^ T d 2 (p A 
— + Th h b h -Q + rPb I-f = 0, 

OX (Jt ox 



* -yp h A h — = 0, M = -E h I^-, Q = KA h G h (—-<p) 



dt 2 



dx 



dx 



(5) 



where E h , G b are the Young's modulus and shear modulus, respectively; ic=7r 2 /l2 is the shear 
correction factor; / is the inertia moment; M , Q are the bending moment, transverse force of 
beam, respectively; (p is the rotation angle of beam due to pure bending; b b ,p b ,A b denote the 
width, density and cross-section area of beam, respectively; the mass ratio y is determined by 
following relation: 

Pb A b +2p p A p + 2p a A i 



r= 

Pb A 

A b =h b b b ,A a =h a b a ,A p =h p b p 

where p a is the density of adhesive layer. 

From Equations (1-5), the following equations can be obtained: 

d 2 u G* h b 3 W d 2 u 

dx h h 2 dx dt 



(6) 



■E b I 



a 4 



w 



■7Pb A \ 



d 2 w ir 2 pld 4 w 



J7 + b b T a ~ ~dx~ + T~d7' ''^ oy KG b ' * J dx 2 dt 2 ^°^dt 2 KG b dt 4 



= 0 



(7) 



(8) 



Assuming the high frequency harmonic voltage loaded on the upper and lower PZT patch are 

V(x, 0 = V(x)e iw and V'(x, 0 = V(x) e' (wr+;r) : 
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then the steady state response of coupled structural system can be expressed as: 

u p (x, t) = Up (x) e im , u s (x, t) = iis(x) Q i(0t , t(x, t) = r(x) e iM , w(x, t) - w(x) e im (Q) 

From Equation (8) and Equation (9), we can obtain: 

- = du p = E b Iha d 4 w h b d 2 w | Ih a yp b o) 2 E b | d 2 w | lh a y 2 _plco 4 h a yp b A^o 2 - 
P dx h^G.dx 4 2 d x 2 h b b b G a ^xG b 'dx 2 h b b b G a xG b h b b b G a (1U) 

where, G a = G a (1 + iC b co) . 

Substituting Equation (10) into Equation (7), the following equation can be obtained: 

d w ,dw , d w , — . 

+ k l — - T + k 2 — — + k 3 w = 0 (11) 



dx dx dx 



where: 



1 2E b Ih a h p h a E p E p KG b E b 
a^PpYPb o) 4 y 2 p b Gjp b co 2 Gjp b (o 2 co 2 yp b A b ofyp p p b h b b b G a p p co 



£ p £ b xG b £ b E p h a E b h p E p h a h p KG b IE b xG b E p 2E p h a E b I 

k ( G*Y 2 pW P v r 2 pW G a yp b A b co 2 | yp b A b co 4 ) 
3 h p h a E p E b KG b E p E b KG b h p h a E p E b I E p E b I 

The solution of Equation (11) can be expressed as: 

w(x) = q e" Ax + c 2 e Ax + c 3 e _/?2X + c 4 e^ x + c 5 e~ Ax + c 6 e Ax , x e [x$ , x 3 " ] (12) 

where, c i •••,€) are the undetermined constant coefficients, ±/? ; - (/ = 1,2,3) are the characteristic 
roots of the following equation: 

^ +k x (3 A +k 2 j3 2 +k 3 =0 (13) 
Substituting the Equation (12), Equation (9) into Equation (10), we can obtain: 

£ V = gl c x e-& x + gl c 2 e& x +g 2 c 3 e-& x +g 2 c 4 e& x +g 3 c 5 e~& x + g 3 c 6 e& x (14) 

where 

= E b I! h pA _K J3 2 + ihrpb® 2 ( E b +V) p2 + { ihr 2 _pW _ KrpbA^ 2 ^ ( i=W ) 

h b b b G a 1 2 1 h b b b G a xG b 1 h b b b G a xG b h b b b G a 
Substituting the Equation (12), Equation (9) into Equation (7), we can obtain: 

-A A -A Pi -A A (15) 

By substituting the Equations (2), (1 1), (15) into Equation (5), respectively, the bending moment M 9 
shear stress Q and rotational angle <p can be obtained. In order to solve the unknown constants 
c. (/ = !,••-, 6), the boundary conditions and the continuity conditions are required, so, the analytical 
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expression of other beam section is also needed. The basic equations of a uniform beam section 
without bonded PZT patches are: 



dM 



d 2 <P_ n dQ 



w 



.dtp 



3 Q + p b I-^ = 0,f--p b A b — = 0,M = -E b I^,Q = KA b G b (—-<p) 

ox at ox at ox ox 



(16) 



2 pllo) A 



From Equation (16), the following relations can be obtained: 

E h I + (—f- + l)//7 b ^ 2 -— - - (p h A h co z - ^—)w = 0 
dx 4 xG h dx z *G h 

The steady state solution of the Equation (17) can be expressed as: 



(17) 



(18) 



w(x) = c 7 e 2 +c 8 e 2 + c 9 e 2 +c 10 e 2 ,ie [0,^ ] 
l i 1 j 1 j 1 j 

— — k 4 x —k 4 x — k 5 x —k$x 

w(x) = c n e 2 +c 12 e 2 +^136 2 +c 14 e 2 ,xe[xi ,x^] 

V -V V 

w(x) = c 15 e 2 +c 16 e 2 +^176 2 +c 18 e 2 ,xe[x 3 ,/] 
where, c,- (/ = 7,. . .,18) are the undetermined constant coefficients; and £ 4 , £ 5 can be expressed as: 

k 4 = {-2k 6 -2^| -4* ? ] , £ 5 = (-2* 6 +2^-4k 7 

KG b E b E b I E b KG b 

Substituting the Equation (18) into Equation (16), we can obtain the bending moment M, shear 
stress Q and rotational angle <p of the beam without bonded PZT patches. In this study, the crack is 

assumed to be a fully open crack, and the depth of crack is relative small. Therefore, the crack on the 
beam can be modeled by a weightless rotational spring [31] as shown in Figure 3. The softness of the 
spring 6 is a function of the beam thickness and the depth of crack [31,32] and can be expressed as: 

0 = 6;r(l-// 2 )/z b (0.6272£ 2 -0.4533£ 3 +4.5948£ 4 -9.9736£ 5 

+ 20.2948£ 6 -33.035 if +47. 1063£ 8 -40.7556£ 9 +19.6£ 10 ) 

where £ = h c /h b ; h c is the depth of crack, and ju is the Poisson's ratio of the beam. 

Figure 3. The equivalent model of the cracked beam. 



(19) 




W////A 
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Considering the continuity of transverse placement, bending moment, shear stress and discontinuity 
of slope at the crack, we can obtain the following equation: 

w + = w -,M + =M,Q + =Q,(p + =(p-+e^- (2Q) 

where, "+" denotes the right side of crack; ' -"denotes the left side of crack. 
2.2. Steady-State Solution 

There are 18 undetermined constant coefficients c t (i = 1,---,18) in the steady analytical solution of 
coupled system. The undetermined coefficients can be determined through boundary condition and 
continuity condition. For the classical end of beam, one has the following equations: 

hinged end 

w = 0,M =0 



clamped end 

w = 0, <p = 0 

free end 

M = 0,g = 0 



(21) 
(22) 
(23) 



For the cantilever beam, the boundary conditions are expressed as: 

w(0) = ^(0) = 0 (24) 
M(O) = 0(O) = O 



(25) 



The boundary conditions of PZT patch can be expressed as: 



£ v (x 2 ) = d 3l E 3 = ^^,£p(x 3 ) = d 3l E 3 = (26) 

/Zp /Zp 

According to the continuity of beam deflection, slope, bending moment and shear stress at the 
borders of the sections at x = x x , x = x 2 , x = x 3 we can obtain the following equations: 

(P{X[ ) = (p{x\) + x - x ~ I <P(*2) = <P(4\ ^>3~) = <p(4l ( 2? ) 

dx 

M(xf) = M(x 1 + ), M(x2) = M(xJ) ? M(x 3 ") = M(x 3 + ), (28) 

Q(xD = Q(4), Q(x2) = Q(4l Q(*D = Q(4) (29) 

w (xf) = w(x 1 + ), w(x2) = w(x2), w(x 3 ) = w(x 3 ) (30) 

The 18 unknown coefficients c t (i = 1 , - - - , 1 8) are obtained from the system of algebraic equations 
determined by the boundary and continuity conditions [Equations (24-30)]. 
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3. Electromechanical Signatures 

Consider a pure extension of PZT patches, and then the corresponding constitutive equations of 
PZT patches can be expressed as [33]: 

£\=°\l E p+ d 3\ E 3 

-T (31) 

D 3 =d 3l a l +€ 33 E 3 

where e x , o x are the strain and stress along x direction, respectively; D 3 , E3 are the electric flux 
density and electric field intensity along height direction, respectively; E = E p (l + jrf) is the complex 
Young's modulus of the PZT material at zero electric field with 77 denoting the mechanical loss factor; 
€33 =€ 33^~ j$) i s complex dielectric constant at zero stress with 5 denoting the dielectric loss 
factor of PZT patch. 

From Equation (31), we can obtain the following equation: 

D 3 = d 3X {s x - d 3X E 3 )E p + e T 33 E 3 = d 3X E p £ x + (4 - 4 x E p )E 3 (32) 

The electric current passing through the upper PZT patch can be determined from the electric 
displacement as: 

Jx 2 



I v= ico l b 0 a r x 3D 3 dxd y 



= icob^d 3X E v [u v {x 3 )-u v {x 2 )] + icoE 3 b^x 3 -x 2 )(£ 33 -d 3X E v ) 

where I p is the electric current passing through PZT patch. 

Electric admittance of the upper PZT patch can be determined as: 

Y = !jl = i(0 K d 3i E v (u p (x 3 )-u p (x 2 )) + (g 3 3 ~ d 3X E v )b a (x 3 - x 2 ) 
V V \ 

= ia) K d 3l E p( U v(X3)- U v(X 2 )) + ia) (^33- d 3l E p) b Jp 

v \ 



(33) 



(34) 



= ico * ?! p £ ^ [(e Ax3 - e ^ )c 2xi - (e^ 3 - ^ )c. ] + i 



CO- 



where 7 is the electric admittance of PZT patch. 

For a piezoelectric system without damage on the PZT patches, the parameters of the PZT patches 
can be regarded as constants. Hence, the change of admittance of PZT patch is only determined by the 
first term on the right side of Equation (34). From Equation (34), we can see that any changes in beam 
structure will lead to a change of the admittance signature of the PZT patch. 

4. Numerical Examples 

In this paper, a cantilever beam with a pair of PZT patches bonded symmetrically onto its top and 
bottom surfaces is studied. The PZT patches are located at X2 = 200 mm from the left side of beam 
(Figure 1). Geometric parameters and material constants of the beam, PZT patches and adhesive layer 
are listed in Tables 1-3, respectively. 
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Table 1. Properties of the beam. 



l h (mm) 


K ( mm ) 


£ b (Pa) G b (Pa) Poisson's ratio /> b (kg/m 3 ) 


Damping ratio 


300 


8 


6.60E10 2.33E10 0.33 2700 


0.01 


Table 2. Properties of the PZT patch. 


/ p ( mm ) 


h p (mm) 


b p (mm) £ p (Pa) /> p (kg/m 3 ) J 31 (m/V) £ 33 (F/m) 


S i] 


10 


0.5 


8 6.1E10 7750 -1.71E-10 1.53E-8 


0.02 0.03 



Table 3. Properties of the adhesive layer. 



l a (mm) 


h & (mm) b & (mm) 


G a (Pa) 


y9 a (kg/m 3 ) 


C b Poisson's ratio 


10 


0.1 8 


1.0E9 


1700 


l.OE-7 0.38 



By considering the inertia terms of PZT patches and bonding layers caused by their motion with the 
beam, the values of y against the beam thickness h h can be calculated and the result is listed in Table 4. 
If the inertia terms are not taken into account, then y = 1 . In order to study the effect of inertia, the 
differences of the admittances of PZT under two scenarios (considering/without considering the inertia 
effect) are compared. To verify the accuracy and the reliability of the proposed analytical model, the 
results are compared to the simulation results obtained by FEA (finite element analysis), which can 
predict well the experimental results [34]. Yang et al. developed a multi-physics simulation method of 
EMI modeling, which can effect direct acquisition of PZT electrical admittance [35]. Yang's method 
was also used in this work. The test specimen is numerically modeled in the ANSYS 10.0 workspace 
as illustrated in Figure 4. The constitutive data, in accordance to PZT-5A, are assigned to the PZT 
patch as given in Table 5. The properties of beam and bonding layer are listed in Tables 1 and 3, 
respectively. The PZT patches are modeled with Solid 5 elements and the bonding layer and beam 
with Solid 45 elements. The sizes of the elements are less than 1.0 mm. An alternating (sinusoidal) 
voltage of 1 V was applied across the PZT patch for excitation. The detailed process can be found in 
the references [34] and [35]. 

Table 4. y vs. h b . 

h h (mm) 2.5 5 10 30 

y 3.321 2.116 1.580 1.193 

Table 5. Constitutive data of the PZT-5 A. 



£(m 2 /N) rf(C/N) e (F/m) 



Sn 


= 5 22 =16.4E-12 


d 3l = 


-171E-12 


£11 


= 1.53E-8 


Sn 


= -5.74E-12 


dn = 


-171E-12 


£ 22 


= 1.53E-8 


Sn 


= 5 23 = -7.2E-12 


J 33 = 


374E-12 


£33 


= 1.50E-8 


S33 


= 18.8E-12 


J42 = 


584E-12 






S44 


= 5 55 = 47.5E-12 


dsi = 


584E-12 






S66 


= 44.3E-12 
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Figure 4. FE model of beam and PZT patches in Ansys. 
Beam sturcture PZT patch 




The results shown in Figures 5-7 suggest that the proposed analytical model provides reasonable 
predictions of the the FEA results, as the major resonance peaks are well predicted. For any beam 
thickness h h , the corresponding frequencies to the peaks of admittance curves considering the inertia 
terms are smaller than those without considering the inertia terms; and the difference between the 
curves is increased with the increasing frequency. When h h =10 mm (Figure 7), the difference between 
curves of y = 1.58 and y = 1.0 is small, however, it is increased with the increasing frequency, and 
the difference of peak frequency is about 0.45 kHz near 96 kHz. From Figure 7, It can be seen that the 
inertia effect cannot be ignored if the frequency is high (>80 kHz). When h\, = 5 mm (Figure 6), the 
difference between curves of y = 2.1 16 and y = 1.0 is more obvious than that of h\> = 10 mm, and the 
difference of peak frequency is about 1.5 kHz near 98 kHz. From Figure 6, it can be seen that the 
influence of inertia can't be ignored when the frequency is higher than 40 kHz. When = 2.5 mm 
(Figure 5), the curves of y =3.321 have a distinct difference with the curve of y = 1.0. The difference 
of peak frequency between y =3.321 and y = 1.0 near 41 kHz is about 0.73 kHz. 



Figure 5. The influence of inertia terms of PZT patches and bonding layers when 
h h = 2.5 mm. (a) 40 kHz-60 kHz; (b) 80 kHz-100 kHz. 



0.00022 
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0.00018 - 
0.00016 - 
0.00014 
0.00012 
0.00010 
0.00008 h 
0.00006 
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0.00002 

o.ooooo f- 

-0.00002 




40 42 44 46 48 50 52 54 56 58 60 
Frequency/ kHz 



0.00007 



0.00000- 




80 82 84 86 88 90 92 94 96 98 
Frequency /kHz 



100 



(a) 



(b) 
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Figure 6. The influence of inertia terms of PZT patches and bonding layers when 
h b = 5 mm. (a) 40 kHz-60 kHz; (b) 80 kHz-100 kHz. 



0.00014 




40 42 44 46 48 50 52 54 56 58 60 
Frequency /kHz 
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0.00025 



0.00020- 



8 0.00015 
c 
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E o.oooioh 

o 
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0.00005 



0.00000 
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90 92 94 96 98 100 



Frequency /kHz 

(b) 



Figure 7. The influence of inertia terms of PZT patches and bonding layers when 
h b = 10 mm. (a) 40 kHz-60 kHz; (b) 80 kHz-100 kHz. 



0.00010 
0.00009 
0.00008 
0.00007- 
0.00006 
0.00005 
0.00004- 
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0.00002 
0.00001 
0.00000 

40 




46 48 50 52 54 
Frequency /kHz 



60 



o 
O 




80 82 84 86 88 90 92 94 96 98 100 
Frequency /kHz 



(a) 



(b) 



As shown in Figures 5-7, when the mass ratio y is l ar g e > the inertia forces of PZT patches and 
adhesive layers can influence the admittance signature greatly, even in the low frequency range, and 
the inertia term must be considered; when the mass ratio y is small, the inertial term also need to be 
considered in the high frequency range. Because the EMI technique uses high-frequency alternating 
current, the proposed EMI model which takes inertial term into account will predict more accurate 
results. 

Consider a beam with a crack located at x = 100. The parameters of the beam are listed in Table 1, 
and the beam thickness hb = 5 mm. the depth of crack is h c , the relative depth £ = h c lh h . To verify the 
reliability of the proposed method, the admittance signatures of PZT patch from Ansys are compared 
with the data obtained by the analytical model. The crack is modeled in the ANSYS 10.0 as a slot with 
a width of 0.1 mm. the size of mesh near the slot is smaller than that of other areas (Figure 8). The 
results are shown in Figure 9. From Figure 9, it can be seen that the proposed method can predict 
admittance signatures of damaged structural as well as the FEA method. In order to study the influence 
of the crack on the admittance signature of the PZT patch, the admittance signatures are calculated 
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when £ = 0.05, £ = 0.1, £ = 0.2 and £ = 0.4, respectively. The corresponding results are shown in 
Figure 10, Figure 11 and Figure 12. The resonant peaks of the admittance signatures shift towards the 
left with increasing crack depth. The influence of the crack increased with increasing frequency. These 
results are in good accord with the experimental phenomena in reference [30] and reference [33]. In 
fact, the change of admittance signature and the decrease of the peek frequencies reflect the decreasing 
local stiffness due to the crack, so the EMI technique can be used to identify the crack damage. From 
the changes of the admittance signature caused by the appearance of damage, both the location and 
quantity of the damage can be identified by using a certain back-calculation algorithm. 

Figure 8. FE model of beam with a crack in Ansys. 



The crack has little influence on the admittance signature in some frequency bands, such as 
45 kHz-47 kHz in Figure 10, 80 kHz-83 kHz in Figure 11 and 191 kHz-193 kHz in Figure 12. The 
frequency band is related to the crack location. The same phenomenon was also found by Youdi [33]. 
The reason is that if the crack is located at the node of a displacement modal, the crack has no 
influence on the modal frequency; otherwise, the modal frequency will reduce with the increasing 
crack depth. Hence, in order to improve the damage identification accuracy, the sweeping frequency 
band needs to contain enough peek frequencies to avoid the error due to the absence of influence of the 
crack around some frequencies. From the signatures of different damage extent, we can calculate the 
RMSD (root mean square deviation) value as follows: 



RMSD 



1=1 ' (35) 



u \2 



LOT) 

i=\ 



where Y is the admittance of PZT, the superscripts d and u denote the signature of the damaged 
structure and undamaged structure, respectively, and n denotes the number of sample points. 

The RMSD values are shown in Figure 10, Figure 11 and Figure 12. From these figures, it can be 
seen that for the same extent of damage, the RMSD value increases with increasing frequencies. In 
other words, a higher frequency is more sensitive to damage than a lower frequency. However, when 
the driven frequency is higher than 200 kHz, the signature can be easily affected by temperature and 
bonding layer, hence, the driving frequencies used in EMI techniques are normally less than 
200 kHz [35]. According to the engineering experience with the EMI technique, the minimum RMSD 
value which can be viewed as a reliable indicator of the existence of structural damage can be set as 
1% [36]. In this work, the crack damage with a depth h c = 0.5 mm = 0.1) can be detected, as the 
RMSD value in the frequency range 180 kHz-200 kHz is greater than 5%. 
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Figure 9. Comparison of admittance signatures against frequency plot between the 
ANSYS simulation and the analytical results, (a) 40 kHz-60 kHz; 
(b) 80kHz-100kHz. 
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Figure 10. The influence of different extent crack damage in 40 kHz-60 kHz. 
(a) the admittance signatures of different damage extent; (b) RMSD value of different 
damage extent. 
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Figure 11. The influence of different extent crack damage in 80 kHz-100 kHz. 
(a) the admittance signatures of different damage extent; (b) RMSD value of different 
damage extent. 
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Figure 12. The influence of different extent crack damage in 180 kHz— 200 kHz. 
(a) the admittance signatures of different damage extent; (b) RMSD value of different 
damage extent. 
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5. Conclusions 



An EMI model of a cracked Timoshenko beam with a pair of PZT patches system has been 
developed by considering the inertial forces of both the PZT patches and the bonding layers. The 
theoretical analysis and numerical tests are focused on the influence of crack and the influence of 
inertial forces because of the beam transverse motion. Through numerical tests results, the following 
conclusions can be drawn: 

(1) The inertial forces of PZT patches and bonding layers produced by the transverse motion of 
beam can be ignored in the low frequency band but should be considered in the high frequency band, 
especially for a thin beam structure. Because the EMI technique employs high frequency, taking the 
inertial forces into account is necessary when monitoring a thin beam structure. 

(2) The admittance signature of the PZT patch can reflect the crack damage very well, especially in 
the high frequency band. In some frequency bands, the crack has little influence on the admittance 
signature, while in other frequency bands the peak frequency of the admittance signature decreases 
with increasing frequency. In order to improve the accuracy of damage identification, the high 
frequency band which contains many peak frequencies should be chosen. 

Based on the proposed EMI model, future work is planned to identify the crack damage 
quantitatively. To quantify the crack damage in plates and shells, further research is needed to 
establish the corresponding EMI models for cracked plates and cracked shells. 
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